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We study the behaviour of the spatial and temporal 
't Hooft loop at zero and finite temperature in the 4D 5(7(2) 
Yang-Mills theory, using a new numerical method. In the de- 
confined phase T > T c , the spatial 't Hooft loop exhibits a 
dual string tension, which vanishes at T c with 3D Ising-like 
critical exponent. 

PACS numbers: 12.38.Gc, 12.38.Aw, 11. 15. Ha 



The AD SU (2) Yang-Mills theory undergoes a transi- 
tion between a cold confined phase and a hot deconfincd 
phase at a critical temperature T c . An order parameter 
widely used to characterize this transition is the Polyakov 
loop. It develops a non- vanishing expectation value in the 
deconfined phase. However, the corresponding operator 
creates a single fundamental static color source, which 
does not belong to the physical Hilbert space of the the- 
ory; it cannot be defined at zero temperature; and it 
is afflicted by ultra-violet divergences in the continuum 
limit [QJ^] . A long time ago it was proposed ||] to consider 
instead the 't Hooft loop operator as an order parameter 
to characterize this transition. It is the purpose of this 
paper to study this order parameter. The 't Hooft loop, 
W(C), is an operator associated with a given closed con- 
tour C, and is defined in the continuum SU (N) theory 
by the following equal-time commutation relations || 

[W{C),W{C')\ = \W(C),W(C')] = (1) 
W\C)W{C')W{C) = e^oc-WiC) (2) 

where W(C) is the Wilson loop associated with the 
closed contour C and ncc' is the linking number of C 
and C". Just like the Wilson loop creates an elementary 
electric flux along C", the 't Hooft loop creates an ele- 
mentary magnetic flux along the closed path C affecting 
any Wilson loop "pierced" by C. In that sense, the two 
types of loop are dual to each other. At zero tempera- 
ture, it has been shown (3[j^] that also the 't Hooft loop 
behaviour is dual to that of the Wilson loop: in the ab- 
sence of massless excitations, an area law behaviour for 
one implies a perimeter law for the other, and vice versa. 
Hence, at T = the 't Hooft loop obeys a perimeter law. 

Several analytical [^|J|] and numerical studies 
have been carried out in order to investigate this issue 
of duality at finite temperature. At T > 0, the Lorentz 
symmetry is broken, so spatial and temporal loops can 



have different behaviours. Because the spatial string ten- 
sion persists also above T c for the Wilson loop, temporal 
't Hooft loops are expected to show a perimeter law in 
both phases; spatial 't Hooft loops are expected to obey 
a perimeter law in the confined phase and an area law - 
defining a dual string tension (strictly speaking it is an 
action density) - in the deconfined phase. 

On the lattice the 't Hooft loop is defined as follows. 
Let us consider the SU (2) lattice gauge theory with the 
usual Wilson plaquette action, 



S(l3) 



/3£(l-iTr(t/ F )), 



(3) 



where the sum extends over all the plaquettes P, and 
Up is the path-ordered product of the links around P. 
Starting from S(/3), one defines the partition function 



Z(I3)= [dU]eM-S(P)) ■ 



(4) 



Let us now switch on "by hand" an elementary magnetic 
flux along a closed contour C defined on the dual lattice. 
To create this magnetic flux, we have to multiply Up by a 
non-trivial element of the center group for the plaquettes 
P dual to a given surface S supported by C H,Q. For 
the SU(2) gauge group, this means flipping the sign of the 
coupling, since the only non-trivial element of the center 
Z2 is —1. Call ViS) the set of plaquettes whose coupling 
is flipped f3 — > — j3. Then the action of the system where 
an elementary flux along a closed contour C has been 
switched on is given, up to an additive constant, by 

Ss(P) = -\p[ £ Tr(U P )- ]T Tr(U P )\ , (5) 



and the partition function is 



Zc{fi) = j W] exp(- 



■Ss(P)) 



(0) 



Zc{f3) does not depend on the particular chosen surface 
S, since different choices are related by a change of inte- 
gration variables. The simplest choice for S is the mini- 
mal surface spanning C. Thus, if C is an R x x R y rect- 
angle in the (x, y) plane (spatial 't Hooft loop), one flips 
the coupling of the (z, t) plaquettes dual to the plaquettes 
belonging to the rectangular area. 
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The 't Hooft loop expectation value is related to the 
free energy cost needed to create the magnetic flux along 
the contour C and is given by 

(W(C)) = Z c (f3)/Z(f3) . (7) 

This expression can be rewritten in the form 

(W(C)) = (exp(-/3 Tr (^))>' (8) 

Pev(S) 

with the average taken with the standard Wilson action. 
In this form, the difficulty of measuring the 't Hooft loop 
becomes clear: the observable is exponentially suppressed 
on typical configurations of the statistical ensemble, and 
gets a significant contribution only from configurations 
having an extremely small statistical weight. Therefore 
the numerical evaluation of (W(C)) represents a difficult 
sampling problem, increasingly so with the loop size. 

Recently, the 't Hooft loop, or special cases of it, 
have been studied numerically on the lattice. In [QJ|] 
the sampling problem was overcome by using a multi- 
histogram method, where one performs several different 
simulations in which the coupling associated with the 
plaqucttcs in V(S) is gradually changed from (3 to — f3. 
In ||, instead, the derivative A/Al3\n(W(C)) has been 
determined, which is a much simpler numerical task. 

In this letter we report on a similar numerical study 
where, by measuring spatial and temporal 't Hooft loops 
at zero and finite temperature, we confirm the role of 
the 't Hooft loop as a dual order parameter for confine- 
ment. We adopt a new numerical method, rewriting the 
ratio Zc{f3)/Z({3) as a product of intermediate ratios, 
each easily measurable. We establish the perimeter law 
behaviour of the 't Hooft loop at zero temperature. We 
measure the free energy of two center monopoles as a 
function of their separation: as expected from duality, 
temporal 't Hooft loops show screening behaviour at all 
temperatures, while spatial 't Hooft loops exhibit a dual 
string tension above T c . Moreover, this dual string ten- 
sion vanishes at T c with a critical exponent v very close to 
that of the 3D Ising model, consistent with universality. 

The numerical technique - We now describe our nu- 
merical method to measure the expectation value of the 
't Hooft operator. Because the relevant contributions 
to Zc{P) and to Z((3) come from regions of the phase 
space with very poor overlap, a direct evaluation of (0), 
(Q) by a single Monte Carlo simulation is not a reliable 
way to compute (W(C)}. It is necessary to consider a 
sequence of intermediate partition functions which in- 
terpolate between Z(j3) and Zc(/3)- The approach used 
in consists of interpolating in the coupling of the 
flipped plaquettes in P(S), from [3 to —0. We inter- 
polate instead in the number of flipped plaquettes, and 
rewrite Zc{(3) / Z((3) as the product of ratios of partition 
functions where the number of plaquettes with flipped 



coupling is progressively reduced. N being the number 
of plaquettes in V{S) 1 we use the identity 

Zc(f3) _ Z N (P) Z N ^{(3) Z x {p) 
Z(p) Z N _ x {p) ' Z N _ 2 {(3) Z (P) [ ' 

where Z k (0), k = 0, . . . , N (Z N = Z c and Z = Z) is 
the partition function of the system where only the first 
k plaquettes in V(S) have flipped coupling. Every ra- 
tio Zk/Zk-i can be now computed efficiently by a single 
Monte Carlo simulation, due to the good overlap of the 
relevant phase space in the two partition functions. Fur- 
thermore, in the practical implementation, it is useful to 
reexpress Z k {(3) / Zk-i(0) as a ratio of expectation values 

Z k (J3) = (exp(-l/3Tr([/ P J)) fc 
Z k _ x (P) (exp(+i/3Tr(C/ P J)) fc 1 > 

The averages (-)k are computed with respect to an action 
where the fcth plaquette Up k of P(S) has zero coupling, 
the first (k — 1) plaquettes have coupling — f3 and the re- 
maining ones coupling The benefits over simply com- 
puting (exp(— (3Tr(Up k ))) using the distribution corre- 
sponding to Zk-i(j3) are two-fold, (i) Using the interme- 
diate weight, (•}&, allows a much better sampling for both 
(exp(-i/3Tr(E/ P J)) fe and <exp(+±/3Tr([/ P J)) fe , with re- 
duced errors, (ii) Since the quantity Zk{(i) / 'Zk-i(P) 
refers to a single plaquette of the lattice, it is very use- 
ful to perform a partial integration after each updating 
sweep of the whole lattice, and measure exp(+f Tr([/p fc )) 
and exp(— ^Ti(Up k )) several times while updating only 
the four links belonging to plaquette P^. Since this pla- 
quette has zero coupling, these four links are decoupled 
from each other. Therefore each link can be updated in- 
dependently iVhit times and the different copies of each 
link can be combined to obtain N^ it measurements of 
exp(±^Tr(t/p f .)). Although these measurements are cor- 
related, the variance reduction is important. 

The statistical error on ratio ( Jio| ) must be evaluated 
with care, since the two averages are computed from the 
same sample of configurations. We use a jackknife analy- 
sis. Then, after each Zk(f3)/Zk-i(/3) has been computed, 
the 't Hooft loop expectation value is evaluated according 
to Eq. ([}]) . The final statistical error is obtained by stan- 
dard error propagation since each Zk(f3)/Zk-i(f3) comes 
from an independent Monte Carlo simulation. 

An advantage of our method over the multihistogram 
technique |t],|| is that the products in (^|) give us infor- 
mation on smaller 't Hooft loops for free; moreover, the 
error analysis is simpler and less delicate than in a mul- 
tihistogram analysis. 

Results - We focus on the free energy F(R) of a pair of 
static center monopoles as a function of their separation 
R. It can be obtained as limp^oc — Ln[(W(R, R t ))]/Rt 
by taking elongated R x R t rectangular loops, in the 
same way as one extracts the static potential between 
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two chromo-electric charges. We take R t as large as pos- 
sible, i.e. equal to the lattice size L. This is analogous 
to measuring the correlation of two Polyakov loops, and 
is the correct approach at finite temperature. Therefore 
we must flip the coupling of R x L plaquettes. We do 
this by scanning first the Rt — , then the R— direction. 
With this ordering, the intermediate partition functions 
Zl, Z2L, ■■, Zr x l in (|^) provide us with the free en- 
ergy at separations 1, 2, .., R respectively. The final ratio 
Zlxl/Zq gives the free energy of a center vortex as com- 
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FIG. 1. Partition function Zk/Zo versus k. Incrementing k 
always increase by 1 the area of the 't Hooft loop, but only 
changes its perimeter by ^2 if k — or 1 mod(lO) respectively. 
Zu is clearly sensitive to changes of perimeter but not of area. 

At zero temperature the 't Hooft loop is expected to 
obey a perimeter law: Zk/Z oc e~ cPk , where i\ is the 
length of the contour Ck ■ For the sequence of flipped pla- 
quettes defined above, i\ = 2L if mod(/c, L) — 0, (2L+2) 
otherwise, unless k < L or k > L(L — 1). should there- 
fore center around two values only. This is exactly what 
appears in Fig.l, on a 10 4 lattice. The free energy is 
clearly insensitive to changes in the 't Hooft loop area. 
This is also confirmed by a direct measurement of Creutz 
ratios x(i?, R), which estimate the force at distance R be- 
tween the two magnetic charges and quickly drop to zero 
as R increases. As with the Wilson loop, the coefficient c 
of the perimeter term is affected by UV divergences. We 
verified that it increases as the lattice spacing is reduced. 

The free energy F(r) can be fitted by a Yukawa form 
^—p — , up to an irrelevant additive constant. However 
the screening mass m is rather large, so that the signal 
quickly dies out. Moreover, consistency of data at dif- 
ferent values of the lattice spacing is only fair. This is 
presumably caused by short-distance lattice distortions 
of the Yukawa potential, which could be explicitly taken 
into account. In any case, Fig. 2 shows the free energy as 
a function of the magnetic charge separation. The curve 
corresponds to a screening mass of 2 GeV. It makes sense 
to match this mass with the lightest gluonic excitation, 
the scalar glueball (~ 1.65 GeV), as attempted in 
but this issue awaits a more extensive numerical study. 
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FIG. 2. Free energy of a static pair of center monopoles as 
a function of their separation, at temperature T < T c . The 
scaling behaviour of the data is mediocre. The curve shows a 
Yukawa potential with screening mass 2 GeV. 

At finite temperature we must distinguish between spa- 
tial and temporal 't Hooft loops. In the case of electric 
charges, the spatial string tension persists above T c , while 
the correlation of time-like Polyakov loops develops a dis- 
connected part showing saturation of the static potential. 
Here the temporal 't Hooft loop should show screening, 
while the spatial 't Hooft loop should obey an area law. It 
is clear that the latter will cost more in free energy, since 
the center vortex created by each flipped (z, t) plaquette 
can only spread over a limited time extent T . 
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FIG. 3. Free energy of a static pair of center monopoles as 
a function of their separation at temperature T > T c . The 
spatial 't Hooft loop shows a dual string tension; the temporal 
't Hooft loop shows a screening mass which increases with T. 

Indeed, this is precisely what we observe. Fig. 3 shows 
the same Yukawa form for a temporal 't Hooft loop at 
all temperatures. The only effect of temperature is to 
increase the screening mass. The fitting coefficient of the 
Yukawa potential is smaller than the short-distance per- 
turbative prediction l/47r(27r/<7) 2 Jl^] , but grows towards 
this value at higher (3. An attempt to fit the data with 

the ansatz Fq+c^— Var, including a linear term, gives 

a dual string tension a consistent with zero. In contrast, 
this linear term is required to obtain an acceptable fit 
above T c for the spatial 't Hooft loop: a dual string ten- 
sion appears. As a practical consequence, the screening 
mass becomes yet harder to determine. It seems little 
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affected by temperature, as Fig. 4 (top) shows, unlike the 
temporal screening mass (bottom) which rises more or 
less linearly with T, much like the glueball excitation 
which it presumably represents Jj"3"| . Precise quantita- 
tive statements about these dependences on T require a 
more accurate numerical investigation. 
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FIG. 4. Screening mass as a function of temperature, both 
in units of T c , as extracted from spatial (top) or temporal 
(bottom) 't Hooft loops. Below T c both coincide. The arrow 
gives the mass of the scalar glueball at T = 0. 

The dual string tension a depends on temperature and 
must vanish at T c . Fig. 5 shows that it does so as a oc 
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The critical exponent v, associated with the 

correlation length £ = <7~ 1//2 , comes out very close to that 
of the 3D Ising model: 0.66(3) vs. ~ 0.63. Indeed this 
should be expected since both models are in the same uni- 
versality class, although the Ising exponent is typically 
extracted from the divergence of £ in the symmetric (con- 
fined) phase. This dual string tension can then be taken 
as order parameter for the restoration of the (magnetic) 
Zjv symmetry, corresponding to cfeconfinement p|Jl4]] . 

Fig. 5 also shows high T, perturbative results. Notice 
that a spatial (x, y) 't Hooft loop of maximal size L x L 
introduces a flipped plaquette in every (z, t) plane, which 
is equivalent to enforcing twisted boundary conditions in 
the (z, t) directions, thus creating a Zn interface with 
tension (crT). This interface tension has been calculated 
to two loops Taking the running coupling g(T) from 
p5[ , one obtains the two curves in Fig. 5 for the leading 
and next order. The numerical data lie in between. 

In summary, by using a dual observable, we have mea- 
sured a dual string (or interface) tension in the deconfmed 
phase. The corresponding correlation length diverges at 
T c with the 3D Ising critical exponent v as expected from 
universality. An extension to ST/ (3) is straightforward. 
In that case, since the transition is first-order, the dual 
string tension will persist all the way to T c . Its value at 



T c has been the object of many studies 16 1. Finally, it 
would be very desirable to measure more precisely the 
screening mass and to clarify its physical origin. 













f 












Z 


— 


' 




/ 

X 
+ 


8 = 2 
8 = 2 
8 = 2 


8 — 
6 - 
5 


_ $ 






o 
o 
□ 
w 


6 = 2 
(3 = 2 
8 = 2 
8 = 2 


45 " 
4 

33 " 
31 


- / 














0. 1 




1 







t = (T - Tc) / Tc 

FIG. 5. Dual string tension, in units of T%, as a function of 
the reduced temperature t. The straight line is a power law 
fit to t < 1. The fitted exponent is 1.32(6), to be compared 
with 2v « 1.26 for the 3D Ising model. The curves show 
the perturbative result, to leading (upper) and next (lower) 
order. 
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